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Abstract. It is shown that in a bilayer excitonic superconductor dissipative losses 
emerge under transmission of the current from the source to the load. These losses are 
proportional to the square of the interlayer tunneling amplitude and independent on 
the value of the input current. The case of quantum Hall bilayer is considered. The 
bilayer may work as a transmission line if the input current exceeds certain critical 
value. The input current higher than critical one induces Josephson vortices in the 
bilayer. The difference of electrochemical potentials is required to feed the load and 
it forces Josephson vortices to move. The state becomes non-stationary that leads to 
dissipation. 

PACS numbers: 73.43. Jn, 74.90.+n 
1. Introduction 

Among the phenomena that demonstrate bilayer electron systems in semiconductor 
heterostructures considerable attention is given to superfluidity of electrons-hole (e- 
h) pairs with components belonging to different layers (see, for instance, pQ). A 
flow of electron-hole pairs in the bilayer is equivalent to two oppositely directed 
electrical currents in the layers. Therefore, the superflow of such pairs is a kind of 
superconductivity. For the first time the effect was considered in Refs. j2j [3] for bilayers 
where one layer is of the electron-type conductivity and the other layer - of the hole- 
type one (electron-hole bilayers). The next important step was the prediction of the 
e-h superfluidity for the systems where the conductivity of both layers is of the same 
type [H [5]. In that case the bilayer should be subjected by a strong perpendicular to 
the layers magnetic field. If the total filling factor of the Landau levels is v — 1 the 
number of electrons in one layer coincides with the number of holes in the other layer 
(the empty states in the lowest Landau level play the role of holes) and the Coulomb 
attraction between electrons and holes results in a formation of bound pairs. Note 
that the description of electron-hole pairing in electron-electron bilayers in a quantized 



Josephson vortex motion as a source for dissipation of superflow of e-h pairs 



2 



magnetic field is close to one developed earlier for the quantum Hall electron-hole bilayers 

la ei ei 

The prediction [H [5] has inspired considerable increase of interest to the study 
of this phenomenon, both theoretically 011011111112111311111151113112! and 
experimentally [HI E21EQ1I2HI22E2E1I23E5]. The results of recent experimental 
investigations of quantum Hall bilayers support the idea on superfluidity of e-h pairs 
in these systems. In particular, in the counterflow experiments a huge increase of 
a longitudinal conductivity was observed [TBI fT9l 120] . Other important observations 
are a large low bias tunnel conductivity [21], that also takes place in bilayers with 
large imbalance of filling factors [26], the Goldstone collective mode [22], the quantized 
Hall drag between the layers [23], the interlayer drag [21] and the interlayer critical 
supercurrent [25] in the Corbino disk geometry. The low temperature properties of 
optically generated indirect excitons in bilayers in zero magnetic field were also studied 
experimentally [27! [28], an d specific features in the photoluminescence spectra accounted 
for the Bose-Einstein condensation of electron-hole pairs have been observed. Recent 
important contribution into this topic is connected with the idea of using two graphene 
layers separated by a dielectric layer [29], [301 ED [32] instead of GaAs heterostructures. 
It is expected [29] that in graphene systems in zero magnetic field superconductivity of 
e-h pairs may take place at rather high temperatures. 

A superfluid state of electron-hole pairs in quantum Hall bilayers can be considered 
as a state with spontaneous interlayer phase coherence between the electrons. The 
coherent phase ip is the phase of the order parameter for the electron-hole pairing, and 
gradient of the phase (p determines the value of the supercurrent in the layer (planar 
supercurrent). 

The bilayer system in the counterflow setup [TH [191 [20] can be used as a double- wire 
circuit that transmits the current from the source to the load. If the transmission is 
provided by superfluid e-h pairs one can expect that such a line will work as a usual 
superconducting transmission line. Nevertheless, genuine superconductivity was not 
achieved in the counterflow experiments. A finite resistance can be accounted for a 
disorder that results in a finite concentration of planar vortices in the bilayer [331 EI]- 
But there is another principal circumstance that may forbid genuine superconductivity 
in such a transmission line. To support nonzero current in the load circuit the difference 
of electrochemical potentials between the layers is required. This difference results in a 
temporal dependence of the phase (p. The state with the spontaneous interlayer phase 
coherence emerges if the layers are situated rather close to each other (at the distance 
less or of other of the magnetic length). Therefore, nonzero interlayer tunneling is 
always present in real physical systems. At nonzero interlayer tunneling amplitude 
the difference of electrochemical potentials results in appearance of an a.c. Josephson 
current and the state becomes non- stationary one. As was shown in [35J dissipative 
losses emerge in the non-stationary state. The situation is different for the loop and for 
the load geometry [36] . The steady state is not possible in the loop geometry Fig. [lb 
(that corresponds to the transmission line [HI HH [20]) but it can be realized in the load 
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Figure 1. The load (a) and the loop (b) setups. 

geometry shown in Fig. []Ja under condition that the difference of the electrochemical 
potential of the layers is tuned to zero. 

The load geometry [36] is an appropriate setup for the observation of the 
superfluidity of e-h pairs but it cannot be used for the transmission of the current 
from the source to the load. In more complicated circuits, e.g. made of a stack of 
bilayers, non-dissipative transmission of the energy is not possible, as well. Indeed, a 
bilayer with zero difference of electrochemical potentials can be shunted at both ends 
without any impact on electrical currents and voltages in external circuits. It means 
that the removal of such a bilayer from the circuit cannot change its working parameters, 
including the power of losses. An example of such shunting is shown in Fig. [2J The only 
setup where superfluid properties of e-h pairs may lead to lowering of dissipation is the 
system with nonzero difference of electrochemical potentials between the layers. While 
non- stationary states are not free from dissipation, the losses can be exponentially small 
and one can say about almost superconducting state. 

In this paper we consider dissipative processes in the bilayer in the loop geometry 
and find the conditions at which the losses are small. The state with a.c. Josephson and 
planar supercurrents that emerges in the bilayer can be described as a moving chain of 
Josephson vortices. Under such a motion two mechanisms of dissipation come into play. 
They are the second viscosity that results in dissipation in the non-stationary regime 
and Joule losses caused by a.c. electrical fields that emerge due to spatial and temporal 
variation of the interlayer voltage. 

In Sec. |2] starting from the microscopic Hamiltonian and using the BCS-like many- 
particle wave function we derive the stationary continuity equation. In Sec. [3] stationary 
vortex states are considered and the problem of the lower critical current is addressed. 
In Sec. H]the state with moving vortices is investigated and the power of losses caused 
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Figure 2. A circuit with a stack of bilayers (a) and an equivalent circuit with shunted 
bilayers (b) 

by such a motion is computed. We find the power of losses is proportional to the square 
of the tunneling amplitude and depends nonlinearly on the resistance of the load. 

2. The model 

Let us consider a bilayer electron system subjected by a strong perpendicular to the 
layers magnetic field B. The filling factors of the layers satisfy the condition v\ + z/ 2 = 1, 
where Vi = 2tt£ Tli* £1X6 the electron densities in the layers, and t = ^hc/eB is 
the magnetic length. Implying the energy gap between the Landau levels hu c (u c is 
the cyclotron frequency) be larger than the Coulomb energy e 2 /e£ (e is the dielectric 
constant) we use the lowest Landau level approximation. In this approximation the 
Hamiltonian of the system has the form 

H = ~t^2 (atx a 2X +H.C.)+— ^ ^ V n,n'{^\pn{^Pni-^-8nn^'^Pn{^)]-{^) 
X ra,n'=l,2 q 

Here a^ x , a n x are the creation and annihilation operators for the electrons in the n-th 
layer, X is the guiding center of the orbit, S is the area of the layer, t is the interlayer 
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tunneling amplitude, p n (q) = J2x a t,x+ q ^/2 a n,x-q v £ 2 /2^ qxX ~ q 1 /4 is the operator for 
the Fourier component of the electron density, and V n y(q) = (27ie 2 /eq) exp(— qd\n — n'\) 
is the Fourier component of the Coulomb potential (d is the distance between the layers). 
In what follows we consider the bilayer of a rectangle shape (S = L x x L y ) with the 
planar currents directed parallel to the x axis. 

The state with the interlayer phase coherence is described by the wave function 

\^) = l[(ua+ x + va+ x )\0), (2) 
x 

where the coefficients u and v satisfy the condition \u\ 2 + \v\ 2 = 1. The physical meaning 
of this function and its applicability for a description of the superffuid state of e — h 
pairs in quantum Hall bilayers was discussed in 0, Ell EH] and in a number of further 
papers. Here we just remind the main points. The function (j2j) can be presented in 
another equivalent form 

l 1 ^) = Y\ ( u ~^ va tx^tx) \ vac )-> (3) 
x 

where = a n is the creation operator for the hole, and the vacuum state \vac) is the 
state in which layer 1 is fully occupied and the layer 2 is empty (v\ = 1 and z/ 2 = 0). 
It was shown [8| that the function ([2]), ([3]) corresponds to the ground state of the 
Hamiltonian ([I]) in the limit d/£ — > 0. One can see that the function ([3]) is the standard 
BCS form for the exciton condensate. The exciton is thought of as an electron in the 
layer 2 bound to a hole in the layer 1. One can show (see, for instance, [37]) that in 
the first quantizied language the function ([2]), ([3]) is reduced to the (111) Halperin wave 
function. At finite djl the state ([2]), ([3]) is not the exact ground state. Nevertheless, 
numerical studies show (see e.g. Ref. [38J ) that at djl < 0.5 the overlap with the exact 
ground state is close to hundred percents. For the system with given filling factors of the 
layers the u — v coefficients in ([2]), ([3]) read as u = ^Jv{ and v = e^^/v?. In the absence 
of the interlayer tunneling the energy of the state ([2]), ([3]) does not depend on the phase 
(p. In case of spatially dependent phase = <f(X) the function ([2]), ([3]) describes the 
state with nonzero counterflow electrical currents in the layers (see [9]). 
The order parameter for the exciton condensate at T = is 

A(X) = |a 2X M*> = u* x v x = A e^ x \ (4) 

where A = \fv\ ( 1 — v\ ) • Due to the absence of the kinetic energy of carriers in the 
Landau level the self-consistence equation for the order parameter at T = is reduced 
to an algebraic equation. At nonzero temperature the equation for the order parameter 
takes a self-consistent form ([HI [39] ) from which the temperature dependence of the 
order parameter and the mean-field critical temperature can be obtained. But due to 
the two-dimensional nature of the bilayer excitons the temperature of transition into the 
superfluid state is not the mean-field critical temperature, but the Berezinskii-Kosterlitz- 
Thouless transition temperature T B kt ~ 7r Ps/2, see [9] (the definition of the superfluid 
stiffness p s is given below). 



Josephson vortex motion as a source for dissipation of superflow of e-h pairs 



6 



If the gradient of the phase is small in comparison with the inverse magnetic length, 
the energy of the system E = (^f\H\^f) can be written in the continuous approximation 



E = E + d 2 r 



dx ) 



t 



2ir£ 2 



'costp — 1) 



where E is the energy of the ground state, t 
A 



Ps 



'Ane£ 



IT 

2 6XP 



erfC 



2tA is the tunneling energy and 
d \ ( d 2 \ d 

1 + ^7 



£V2 



(5) 



(6) 



is the energy parameter called the superfluid stiffness. The superffuid stiffness can also 
be presented in a more familiar form p s = h 2 n s /mB, where n s is the superfluid density 
of the pairs, and is the magnetic mass of the pair (see, for instance, [33]). One 
can see from (jSJ) that at nonzero tunneling amplitude the ground state corresponds to 
the phase <p = 0. But it does not mean the fixation of the phase and the absence of 
electrical currents. It means that the planar current should be accompanied with the 
interlayer current. 

Varying the energy ([5J) with respect to the phase and equating the result to zero 
we obtain the equation 



ep s d 2 p> e 



t 



(8) 



n dx 2 1 h2,p s]nip = - (7) 

The first term in Eq. ([7]) is the 2D divergence of the density of the planar supercurrent 

ep s dip 

ji = -32 = — rir- 

a dx 

The expression for the planar current can be obtained from the gauge invariance 
arguments (see [T5| BO]). 

One can see the Eq. ([7]) is the stationary continuity equation 

djl + T - 
-; r i\^2 — 

dx 

where 

t 



0. 



(9) 



sm<p 



(10) 



h2Tt£ 2 

is the density of the interlayer supercurrent flowing from the layer 1 to the layer 2. 

The coherent electrical current between the layers corresponds to nonzero <p, in 
a close analogy with the Josephson effect between two bulk superconductors. The 
difference is that <p is the phase of a given condensate, but not the phase difference for 
two condensates. 

Eq. (J7J) can be presented in the form 

d^ = P sm ^' 

where A = £\/2np s /t is the Josephson length. As follows from Eq. pip , the gradient 
of the phase is of order of A -1 . Therefore, the continuity approximation (jSJ) requires 
A ^> £. This inequality is fulfilled if the tunneling amplitude t is much smaller than the 
Coulomb energy e 2 /e£ and the filling factor v is not very close to zero. 
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3. Stationary vortex state and lower input critical current 

The problem of the lower critical current was already discussed in a number of papers 
[T2| [T5l W\\ . Here we revisit this question again with the aim to give a more accurate 
definition of the lower critical current and to clarify some discrepancies in [321 I4T] - 
Let us consider the following situation. A fixed input current from a source is 
entered into one layer in a given (source) end of the system. There is no difference 
between the electrochemical potentials of the layers and the stationary state is realized. 
The current withdrawn from the adjacent layer at the same end is equal to the input 
current. The values of the output and input currents at the opposite (load) end are also 
equal to each other, but, in general, they may differ from currents at the source end. 
For instance, such a situation corresponds to the load geometry |36j (then all input and 
output currents are equal to one another). In the loop geometry the stationary state 
can be realized if the load has zero resistance (superconducting load), or if an additional 
source in the load circuit provides zero difference of electrochemical potentials between 
the layers. 

Depending on the values of the input currents two qualitative different stationary 
current patterns can emerge: either the planar and interlayer currents are nonzero in the 
whole bilayer, or the currents decrease exponentially with the distance from the ends 
and there is no currents in the internal part of the bilayer. We define the lower critical 
input current as the current at which the switching between these two regimes occurs. 
If the input currents at both ends are lower than the critical one the second pattern is 
realized and the interior part of the bilayer is not involved into the transmission of the 
current (the cutting of the bilayer does not change currents in outer circuits). If one or 
both input currents exceed j c the first pattern emerges. 

Let us switch to the quantitative analysis. Eq. (fTTT) coincides in form with the 
equation of motion for a nonlinear pendulum 



where ujq is the resonant frequency of the pendulum, and (p is the angle coordinate 
counted from the unstable equilibrium point. The formal coincidence of Eq. (Till) and 
Eq. ffT2]) allows us to describe the stationary states in the bilayer basing on known 
behavior of a nonlinear pendulum. 

Depending on the energy, the pendulum oscillates, completes a full revolution in 
infinite time, or rotates. The maximum angular velocity oo m of the pendulum (velocity 
at (p = 7r) is proportional to the square root of the energy and it increases if one switches 
from the oscillating to the rotating regime. For the full revolution in infinite time the 
maximum angular velocity u m = 2ujq. 

The time dependent angular velocity in the pendulum problem corresponds to the 
space dependent planar supercurrent in the bilayer problem. The counterpart of the full 
revolution regime is a state with one Josephson vortex. In this state the planar and the 





(12) 



dt 2 
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interlayer supercurrents are given by the expressions 




) 



(13) 



where x$ is the vortex center, and 

is the maximum value of planar supercurrent in the single vortex state (it is reached at 
the center of the vortex). Note that Eq. (j!4p can be obtained directly from Eq. (181 
under accounting the correspondences (d(p/dx) max <=> u m = 2cj and loq •<=>■ 1/A. 

The counterpart of the rotation regime is the state with many equally distanced 
Josephson vortices with the same sign of vorticity. The supercurrents in this state read 
as 



The analog of the oscillating motion of the pendulum is the multivortex state with 
vortices of alternating vorticity: 

ji = ~h = Jcv^cn IT' 7 ?) > 



In Eqs. ( Il5j) . ffTBl) dn(u,rj), cn(u,rj) and sn(u,7]) are the Jacobi elliptic functions and 
the parameter r\ belongs to the interval (0, 1). 

The boundary condition for the bilayer problem is the condition on the gradients of 
the phase. This condition alone does not determine an unique current state. The system 
chooses the state that at given input currents has the lowest energy. Direct calculation 
shows that the energies of the multivortex state ffl5|) and ffTBT) are both higher than 
the energy of the single vortex state ( ITBl) . It is obvious result since the energy of any 
multivortex state is proportional to the number of vortices. 

Here we should emphasize the difference between the pendulum problem and the 
bilayer problem. The energy of the pendulum is the integral of motion of Eq. ( 1121) . 
while the bilayer energy is not the first integral of Eq. fTTTl) . Therefore, the analogy 
between these two problems fails if one analyzes the energies of different states. 

One can see that the planar supercurrent may exceed j c only in the state ([15]) . 
Therefore, if the input current j m is higher than j c the multivortex state (1T51) is realized. 
The parameter rj is determined by the conditional minimum of energy and it is equal 
to r\ — (jc/jin) 2 - The density of the Josephson vortices increases under increase of the 
input current. 

If the input current is smaller than j c the state with one incomplete vortex centered 
outside the bilayer satisfies the boundary condition and corresponds to the minimum of 




(15) 




(16) 
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energy. Thus, the state realized at j in < j c is the state with zero supercurrents inside 
the bilayer. 

The vortex state with alternating vorticities ( |T6l) is not realized in the bilayers. It 
may satisfy the boundary conditions only at j in < j c , but its energy is larger than the 
energy of the state with one vortex. It contradicts the conclusion of Ref. [41J, but in 
Ref. [H] the first integral of Eq. (ITTj) was incorrectly identified with the energy. 

Thus, the quantity (fli]) is just the lower input critical current defined above. 

At j in ^> j c the overlapping between the vortices is large and the current pattern 
is approximated by a sum of large d.c part and small harmonic a.c. part: 

(p m 7t — kx — A v sin kx, 

Ji = -32 ~ 1 3 ™ A (1 + A v cos kx) , 

e t 

J i-2~T^sinfcr, (17) 



h2n£ 2 



where 



^ ^in 



Here X(?7) is the complete elliptic integral of the first kind. 

To conclude this section we note that at small input current (lower than j c ) the 
tunneling provides shunting of the bilayer at both ends. It means that the load geometry 
experiments proposed in [36] for the observation of e-h superfluidity should be done at 
large currents. 



4. Vortex motion and dissipation 

As was already discussed in the introduction, the stationary state cannot be realized 
in the bilayer transmission line with nonzero load resistance. But even in case of 
non- stationary currents inside a bilayer the currents in external circuits may remain 
stationary one. The latter occurs if the contacts between the layers and the external 
wires average the input and output currents. 

In such a situation the term "lower critical current" slightly changes its meaning. 
The average current depends on the density of vortices and can be smaller than j c . But 
the input current higher than critical one is required at the beginning of the process to 
induce a chain of moving vortices. 

To study the non- stationary vortex state we will use the set of equations for the 
phase of the order parameter cp and for the local difference of the electrochemical 
potentials between the layers eV, where V is the local interlayer voltage. The first 
equation comes from the non-stationary continuity equation 

dn (9?? djf . . 
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Here n is the excess local density of electron-hole pairs. It is connected with the local 
voltage by the capacitor equation 



f = C, (19) 



where 



C = TTJ f 7 xT (2°) 

4^1-^[l-exp(^)erfc(^)_ 

is the effective capacity of the bilayer system per unit area. The effective capacity 
takes into account the exchange interaction between the layers and differs from the 
classical capacity by the additional factor (the second factor in Eq. (1201 ). It can be 
obtained from the dependence of the energy on the filling factor imbalance (see, e.g. 
[39J). The superscripts "s" and "n" in Eq. ( fl8l) stand for the planar supercurrent 
and for the planar current of quasiparticles. The current of qusiparticles is taken into 
account in Eq. (fl8|) because it is induced by a.c. planar electrical fields that emerge 
in the non- stationary state. These fields equal to E\ = —E2 = —(l/2)dV/dx and the 
quasiparticle current can be presented in the form j[ n ^ = —(o~ n /2)dV/dx, where a n is 
the quasiparticle conductivity. We do not take into account the normal component of 
the interlayer current because the normal tunneling is suppressed due to discreteness of 
the Landau levels. 

The second equation on (p and V can be derived from the equations of superfluid 
hydrodynamics [12], namely, from the equation for the superfluid velocity with 
dissipative terms. The linearized version of this equation reads as 

= -V/i m + V [CsV ■ n s (v s - v n ) + C 4 V ■ v n ] , (21) 

where fi m is the chemical potential per unit mass [i m = fi/m (m is the mass of the Bose 
particle), n s is the superfluid density, v n is the velocity of the normal component, £3, 
C4 are the second viscosity parameters. To apply Eq. (T2TT) for the description of the 
bilayer electron-hole superfluidity one should replace the divergence of the superfluid 
flow in (T2~T1) with the divergence of the planar supercurrent plus the tunnel supercurrent: 
V(n s v s ) =>• (l/e)(djl/dx + and substitute the difference of the electrochemical 

potentials instead of the chemical potential n =^> eV. Using Eqs. ( fTOl . (JE|), and the 
expression for the superfluid velocity v s = —KVip/m (the negative sign is because the 
phase ip describes the coherence between the electrons), we obtain 



(22) 



The quantity 0(v n ) in Eq. (122]) stands for the term caused by the motion of the normal 
component. Below we will neglect the 0(t> n ) term. One can show that this term yields 
the correction to the power of losses quadratic in the dissipative coefficients, while the 
leading term is linear in these coefficients. Since we are interested in the vortex solution, 
the gradients in ( f22l can be omitted. 
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Thus, we obtain the following set of equations for if and V 
dV e d\ e t . a D d 2 V 

c ~m = h Ps d^ ~ 8111 * + y^' (23) 



h^ = eV + a, (stt^a^ - f sin v ) . (24) 



The coefficient = (3m/2ii£ 2 h in Eq. ( J241) is the second superfluid viscosity expressed 
in dimensionless units (with m = m#, the magnetic mass of the e-h pair). Eqs. (I23p . ( |2~4"1) 
contain two terms responsible for the dissipation. The term proportional to <r n describes 
the Joule losses caused by a.c. quasiparticle currents, and the term proportional a v - the 
losses caused by the second viscosity. The equations with the same form of dissipative 
terms were used in [13]. Eq. (T53]) was also obtained in Ref. [H] for the electron- hole 
bilayers in zero magnetic field in the "dirty" limit. 

The channels of dissipation caused by the second viscosity and quasiparticle currents 
are present in any superfluid system, but they do not result automatically in dissipative 
losses. One can see that in the stationary regime the dissipative terms in Eqs. (I23l) . ( 1241) 
equal to zero. The losses appear under presence of additional factors. For the bilayer 
transmission line the main factor is non-stationarity of the interlayer and the planar 
supercurrents. 

As was shown in the previous section, the input current may induce Josephson 
vortices in the bilayer system. At nonzero difference of the electrochemical potentials 
the vortices begin to move from the source to the load. Neglecting for a moment the 
dissipative terms we obtain from Eqs. (123]) . fl24l) the following equation for the phase 

W~d^ = sm ^ (25) 

where the dimensionless variables £ = x/X and r = t/r (r = {M/e)\j2nC]t) are used. 
Eq. fl25l) has the solution <£>(£, r) = — f3r), where (3^1 and <p(u) satisfies the 
equation 

9^ = T^ 81110 (26) 

that up to the notations coincides with Eq. (fTTj) . Using the results of Sec. [3] we obtain 
the following expression for the planar current and for the interlayer voltage 

2e 

jl(x,t) = —p s kdn(kx-ut,r)), (27) 
V(x, t) = — dn (kx - wt, rj) , (28) 



where 



k= , p =, u= . ? (29) 



Eq. ( 1271) describes the moving vortex lattice. The velocity of motion is 
u e 2 (V) e 2 o,RL„ 



Josephson vortex motion as a source for dissipation of superflow of e-h pairs 12 



Here and below (. . .) means the time average. In Eq. (1301) R is the resistance of the 
load (we imply the Ohm law (V) = RL y (jf) for the load circuit). One can see that the 
vortex velocity is proportional to the load resistance and does not depend on the input 
current. The parameter (3 = u v t /X is also proportional to the load resistance. 

Let us switch to the case of small but nonzero dissipation. The power of losses can 
be obtained from the common expression for the Joule losses 

AP = L y l£ dx I X ^V + (jt + jT) (-^) \ ■ (31) 

Integrating (I3~TT) by parts and taking into account the continuity equation ( fT8l) we obtain 
the obvious result that the power of losses is the input power minus the output power: 

AP=(j in - Jout )V L y . (32) 

Here j out is the output planar current, and Vo is the voltage in the output circuit (that 
coincides with the voltage applied to the bilayer system at the input end). In deriving 
(1321) we take into account that (V(dV/ dt)) = for any function V(x,t) periodic in t. 

The difference between the input and the output average currents emerges if the 
average value of the interlayer current / = (h-,2) differs from zero. In what follows 
we consider the case where the power of losses is much smaller than the input power. 
In this case one can neglect the dependence of I on x and approximate the difference 
Jin — jout as jin — j ou t = L X I. Then the power of losses reads as 

AP = IV S. (33) 

The quantity I can be found from the solutions of Eqs. ( l23l) . ( j2~4l) . Here we consider 
the case of large input current j m 3> j c - In this case one can seek for a solution of Eqs. 
(O, (|2lD m the following form 

(p(x, t) = ujt — kx + A v sm(ut — kx) + cos(ut — kx), (34) 
V{x, t) = 1^(1 + A v sin(cut - kx) + By cos(cut - kx)) , (35) 

where u = eVo/h, and the quantity k is connected with the planar supercurrent by the 
relation (jf) = ep s k/h. Substituting Eqs. ($M§, (US into Eqs. ([23]) ([24]) and neglecting 
the terms quadratic in t we obtain the coefficients A^, B^,, Ay, By in the leading order 
in t: 

I (1 - (3 2 + 2nek 2 a v d){\ + 27r£ 2 /e 2 ^a) + (a v ^ + a^) 

A * = ~27tW P s D " ~ ~' (36) 

I af p {l + 2M 2 a v ~a) + a^ 

^"2vr£ 2 P Ps D ' (37) 
t ^fVc ^^V~ 

Av = ~2^Wi s D ~' (38) 

t 1 + 27i£ 2 k 2 a v & - j3 2 

Bv = ~2^i 8 D ' (39) 
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where 



D = (1 - (3 2 + 2Ttt 2 k 2 a v ~a) 2 + (a v ^- + a^j , (40) 



and the following notations are used 



,2 _ CV , r e T f Ps x 0- n h 



ftF"' "' c =2^ V ' = e' S = W' (41) 
Note the equivalence of (3 in (14~TI) and the definition of (3 given above. 

One can see that at t — all the coefficients A^, B^, Ay, By are equal to zero that 
reflects the absence of vortices in the bilayers with zero interlayer tunneling. 

Using Eqs. ( jTUi) and ( |34|) one finds the average value of the interlayer current 

— Gt , . Gt . . , , 

(sin^TT^V (42) 



Since B^ is proportional to t the average interlayer current is proportional to the 
square of the matrix element of the interlayer tunneling. The higher order corrections 
to B v yield the contributions oc t 3 into /. Therefore, for obtaining the power of losses 
in the leading order it is enough to use the linear in t approximation for the phase. 

Using Eqs. (1421) and (1371) we obtain the following expression for the power of losses 

AP t 2 ( B 

■ ' ' (43) 



S Air£ 2 h \ R 



to 



where 



(44) 



Here j = ep s /hi (j ^> j c ) is the current of order of the higher critical current [U] (the 
current above which superfluid state is destroyed). The quantity Ro = h/eL y ^p s C is 
the resonant load resistance. One can see that dissipative losses increases considerable 
at R approaches to R (the condition R = R is equivalent to (3 — 1). We emphasize 
that the approximation solution ( |34|) and the answer f j44|) are not valid in a resonant 
case and at R ~ R they describe the situation only qualitatively. 

Under obtaining Eq. (|33l) we take into account the relation (jl) = j = fp s & and the 
Ohm's law Vq ~ R\L y j in the load circuit. In Eqs. (l4"4"|) the terms in nominator quadratic 
in dissipative coefficients are omitted. 

Thus we conclude that the power of losses is proportional to the square of the 
tunneling amplitude and depends nonlinearly on the load resistance. At small load 
resistance R -C Ro the dissipation is connected in the main part with the conductivity 
of quasiparticles. It is proportional to the square of the load resistance 
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and vanishes at R = . At large load resistance R 3> Ro the main contribution the 
into the power of losses comes from the second viscosity. The power of losses at large 
R approaches to the constant value 

AP t 2 



S ~ 47rm a<f - (46) 
The resonant resistance R does not depend on the tunneling amplitude but it can be 
tuned by a change of d/i (it is the increase function of that parameter). At djl — 1 the 
resonant resistance is approximated as 
h 2 I 
e 2 A L y 

At low load resistance the "bottle neck" of the load circuit are the contacts where 
the interlayer phase coherence is broken. We evaluate from (|4T|) that at longitudinal 
resistivity of the layer p xx = IkQ the resistance of the arm R ~ Ro correspond to the 
length of contacts ~ 10 2 £ ~ 1/im. Since it is rather small length we conclude that the 
experimental situation corresponds most probably to the case R 3> Ro- 

In conclusion we note that in the non-resonant regime the power of losses does not 
depend on the input current. Since the input power is proportional to the input current, 
the efficiency factor increases under increase of the input current. 



5. Conclusion 



We have shown that in the bilayer system with superfluid electron-hole pairs the 
dissipation appears under the transmission of the current from the source to the load. 
The effect is connected with that the electrical currents inside the bilayer becomes non- 
stationary at nonzero difference of the electrochemical potentials between the layers. 

The non-stationary state can be interpreted as the state with moving Josephson 
vortices. It it well known that in type II superconductors the motion of quantum vortices 
results in dissipation, but the latter can be eliminated by pinning of the vortices. In 
this connection one can think that the pinning may also suppress the dissipation in the 
bilayers. But it is not true. The difference between type II superconductors and the 
bilayers is the following. In type II superconductors nonzero voltage along the direction 
of the supercurrent emerges due to the vortex motion. In the bilayers nonzero interlayer 
voltage is required to support electrical current in the load circuit and this voltage causes 
the motion of Josephson vortices. The state is stationary only at zero interlayer voltage 
at which there is no current in the load circuit. 

One can also arrive at this conclusion in another way. The energy of a single vortex 
reads as 




According to Eq. (T48|) the pinning may occur due to spatial variation of the tunneling 
amplitude. In the latter case an approximate solution for the phase can be found by 
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the same way as was done in previous section. The only difference that the coefficients 
A v , B v , in Eq. (jMj) should be replaced with x-dependent quantities. The solution 
Eq. ( jMj) with spatially dependent amplitudes A v and B v also corresponds to a non- 
stationary state. It differs from the state described in the previous section by that the 
phase remains non-stationary in any reference frame, including the frame in which the 
vortices are at rest. One can show that in this case the dissipative losses are also nonzero 
and proportional to dissipative coefficients a v and a. 

Thus, nonzero interlayer tunneling results in that the electron-hole pairs in the 
bilayer cannot transmit electrical energy without dissipation. Nevertheless, since the 
power of losses is proportional to the square of the tunneling amplitude and this 
amplitude depends exponentially on the interlayer distance one can expect that it is 

possible to create bilayer systems with negligible small dissipation. 
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